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ABSTRACT 

In this paper we present a new method for deriving UV-through-IR extinction curves, based on 
the use of stellar atmosphere models to provide estimates of the intrinsic (i.e., unreddened) stellar 
spectral energy distributions (SEDs), rather than unreddened (or lightly reddened) standard 
stars. We show that this "extinction- without-standards" technique greatly increases the accuracy 
of the derived extinction curves and allows realistic estimations of the uncertainties. An additional 
benefit of the technique is that it simultaneously determines the fundamental properties of the 
reddened stars themselves, making the procedure valuable for both stellar and interstellar studies. 
Given the physical limitations of the models we currently employ, the technique is limited to main 
sequence and mildly-evolved B stars. However, in principle, it can be adapted to any class of 
star for which accurate model SEDs are available and for which the signatures of interstellar 
reddening can be distinguished from those of the stellar parameters. We demonstrate how the 
extinction-without-standards curves make it possible to: 1) study the uniformity of curves in 
localized spatial regions with unprecedented precision; 2) determine the relationships between 
different aspects of curve morphology; 3) produce high quality extinction curves from low color 
excess sightlines; and 4) derive reliable extinction curves for mid-late B stars, thereby increasing 
spatial coverage and allowing the study of extinction in open clusters and associations dominated 
by such stars. The application of this technique to the available database of UV-through-IR 
SEDs, and to future observations, will provide valuable constraints on the nature of interstellar 
grains and on the processes which modify them, and will enhance our ability to remove the 
multi-wavelength effects of extinction from astronomical energy distributions. 

Subject headings: ISM:dust,extinction — stars:atmospheres, abundances — methods:data analysis 



1. INTRODUCTION 

A detailed determination of the wavelength de- 
pendence of interstellar extinction, i.e., the ab- 
sorption and scattering of light by dust grains, is 
important for two very different reasons. First, 
since it is a product of the optical properties of 
dust grains, extinction provides critical diagnostic 
information about interstellar grain populations 
(including size distribution, grain structure, and 
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composition), providing guidance for interstellar 
grain models. Second, the accuracy to which the 
intrinsic spectral energy distributions (SEDs) of 
most astronomical objects can be determined de- 
pends on how well the effects of extinction can be 
removed from observations. In both fun- 
damental issue is how accurately the wavelength 
dependence of extinction can be measured. Conse- 
quently, it is essential to have a firm grasp of how 
measurement errors can affect the determination 
of this wavelength dependence. 

This paper is the culmination of a series of 
"techniques" papers published over the past five 
years (Fitzpatrick & Massa 1999; Massa & Fitz- 
patrick 2000; and Fitzpatrick & Massa 2005; 



hereafter FM99, MFOO, and FM05, respectively) 
whose aim has been to develop a technique to 
simultaneously determine the wavelength depen- 
dence of extinction (to higher accuracy than pre- 
viously possible) and the physical properties of 
a reddened star. It represents a continuation of 
our earlier series on the properties of UV extinc- 
tion (Fitzpatrick & Massa 1986; 1988; and 1990, 
hereafter FM90), and provides a detailed descrip- 
tion of a new technique for deriving interstellar 
extinction curves which does not rely on obser- 
vations of standard stars, virtually eliminates the 
effects of "mismatch" error, and yields an accurate 
assessment of the uncertainties. This "extinction- 
without-standards" technique opens the door to 
a new class of extinction studies, including re- 
gions heretofore inaccessible. For example, er- 
rors in the traditional "pair method" approach 
to extinction strongly limit our ability to study 
extinction in two important regimes. The first is 
low-.E(-B — V) sightlines, where one might hope 
to relate extinction properties to the environmen- 
tal properties of specific physical regions. The 
second is extinction derived from mid- to late B 
stars. These stars are plentiful, and often consti- 
tute the bulk of the stars available for extinction 
measurements in intrinsically interesting regions, 
such as the Pleiades. We will demonstrate how the 
extinction-without-standards approach overcomes 
both of these problems and present examples of 
each. Some first results from this program were il- 
lustrated by Fitzpatrick (2004, hereafter F04). In 
addition to the extinction results, we will demon- 
strate that our approach simultaneously provides 
accurate stellar parameters which are also astro- 
physically interesting. 

In §2, we provide a broad overview of the ba- 
sic problem of measuring an extinction curve, and 
compare the merits of curves derived by the pair 
method and curves derived by using stellar atmo- 
sphere models. In §3, we describe our new model- 
based technique in detail, and list the basic ingre- 
dients needed to determine an extinction curve us- 
ing this approach. In §4, we describe the data used 
in the current study. In §5 we provide a number 
of sample extinction curves derived using model 
atmospheres. We also demonstrate the high pre- 
cision of the new curves and the reliability of the 
error analysis employed. Finally in §6 we describe 
some to the scientific advantages of this new tech- 



nique and our plans to exploit them. 

2. MEASURING EXTINCTION 

To understand how interstellar extinction is 
measured and to appreciate how different mea- 
surement techniques can affect the outcome, we 
begin with the intrinsic elements of an uncali- 
brated observation of a spectral energy distribu- 
tion (SED) of a reddened star obtained at the 
earth, f\. This can be expressed as 

h = F x r x 9 2 R W- 0AA > (1) 

where F\ is the intrinsic surface flux of the star 
at wavelength A, r A is the response function of the 
instrument, 9r = (R/d) 2 is the angular radius of 
the star (where d is the stellar distance and R is 
the stellar radius), and A\ is the absolute attenua- 
tion of the stellar flux by intervening dust (i.e., the 
total extinction) at A. Alternatively, the observed 
SED can be expressed in terms of magnitudes ra\ 

by 

m A = -2.51ogF A ^ + A A + C A or (2) 
m A = M A + 51ogd-5 + ,4 A + C A (3) 

where C\ — — 2.5 logr A is a term which transforms 
between the observed magnitude system and abso- 
lute flux units, and M\ is the traditional definition 
of the absolute magnitude of the star at A. 

The difficulty in measuring the total extinction 
A\ can be seen by rearranging these equations to 
solve for the extinction term. Equation 1 yields 

A A = 2.51og(r A ^) , (4) 

/A 

while Equation 3 yields 

Ax = m A - M\ - 5 logd + 5 + C A . (5) 

In either case, a true measurement of A\ would 
require calibrated SED observations, knowledge of 
the intrinsic SED of the star, and measurements 
of both the stellar distance and radius, or their 
ratio Or. Unfortunately there are no early- type 
stars for which both of these latter quantities are 
known to sufficient accuracy to allow a meaningful 
measurement of A\. As a result, indirect meth- 
ods must be employed, and A\ is always a derived 
quantity, subject to assumptions. 
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In virtually all extinction studies, the actual 
measured quantity is a "color excess" which de- 
scribes the extinction at a wavelength A relative to 
that at a fiducial wavelength. The traditional ap- 
proach is to adopt the V band as the fiducial, since 
V magnitudes are widely available, accurately cal- 
ibrated, and typically of high quality. This color 
excess can be expressed as either 

E(X-V)=A x -A v = 2.5log(^^^f) , 

(6) 

as based on Equation 4, or 

E(X-V) = A x -A v 

= (m A - m v ) - (M\ - M v ) + (C x - C v ) 
= m{\-V)-M(\-V) + C{\-V) (7) 

as based on equation 5. Thus, the determina- 
tion of the color excess requires only the measure- 
ment of the observed SED and a knowledge of the 
shape of the intrinsic SED. There are two basic 
approaches to determining color excesses, based 
on the use of either unreddened stars or stellar at- 
mosphere models to represent the intrinsic SEDs 
of reddened stars. These two techniques, and the 
issue of the normalization of extinction curves, are 
discussed in the three subsections to follow. 

2.1. The Pair Method 

The first approach is the "pair method." A 
pair method curve is constructed by comparing 
the fluxes of a reddened star and an (ideally) iden- 
tical unreddened "standard star." Essentially, the 
absolute magnitudes in Eq. 7 are replaced by the 
observed magnitudes of the standard star and the 
curve is usually expressed in the form 

M\-V\ = E ^ X - V ) = ™(.\-V)-m(\-V) 
{ '~E{B-V)~ (B — V) — (B — V)o 

(8) 

where the color excesses E(X — V) arc normal- 
ized by E{B — V) and the subscripted quantities 
refer to unreddened indices for the standard star. 
(The issue of normalization will be discussed be- 
low.) When the two stars are observed using 
the same instrument, this method has the advan- 
tage that calibration terms cancel, eliminating any 
dependence on the absolute flux calibration, r\ 
or C\. There are, however, two disadvantages 
to this technique. The first is that the grid of 
unreddened standard stars is necessarily limited, 



so that some mismatch in the SEDs of the red- 
dened and unreddened stars, termed "mismatch 
error", is inevitable. The second is that there are 
very few truly unreddened early-type stars, so usu- 
ally the "unreddened" standard must be corrected 
for some small amount of extinction whose exact 
magnitude and wavelength dependence are uncer- 
tain. This creates an error that can propagate into 
the resulting extinction curves. 

Massa, Savage, & Fitzpatrick (1983) presented 
a detailed study of the uncertainties affecting pair 
method extinction curves and showed that mis- 
match effects dominate the error budget. If one 
uses a single unreddened spectral standard for 
each spectral class, and assumes that all spectral 
classifications are perfect, then as a result of spec- 
tral binning, mismatch errors will be equal to or 
less than half a spectral class. Figure 1 shows how 
such mismatches can affect extinction curves de- 
rived from main sequence B stars. In each of the 
four groups of curves in Figure 1, the true shape of 
an extinction curve affecting a Bl, B2, B5, or B9 
star is indicated by the solid curve. The extinction 
curves which would be derived via the pair method 
in the presence of a ±| spectral class mismatch er- 
ror are shown by the dash-dot curves for the case 
of E(B — V) = 0.15 and by the dashed curves for 
E(B - V) = 0.30. (The derived curves fall below 
the true curve in the UV and above the true curve 
in the IR when the standard star is cooler than the 
reddened star, and vice versa.) In practice, spec- 
tral classifications are not perfect and the available 
unreddened standard stars do not necessarily lie in 
the middle of the range of properties within a sin- 
gle spectral class. Both these effects exacerbate 
the spectral mismatch problem and thus the un- 
certainties shown in Figure 1 arc likely closer to 
typical errors, rather than extremes. 

Although it may appear from Figure 1 that 
a discontinuity at the Balmer jump at A -1 ~ 
2.7^m~ 1 would provide an obvious indicator of 
the presence of spectral mismatch in an extinc- 
tion curve, this is almost never practical. The 
spectrophotometric data available for construct- 
ing extinction curves are generally limited to UV 
wavelengths (A" 1 > 3.3/im _1 ), which are high- 
lighted in Figure 1. Typically, the only data avail- 
able in the optical and near-UV are photometric 
indices which straddle the Balmer jump, such as 
the Johnson U — B color, and these cannot be used 
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to distinguish the effects of mismatch from intrin- 
sic curve shape. 

Mismatch error clearly can have a profound ef- 
fect on the shapes of curves derived from stars with 
low color excesses, particularly in the mid-to-late 
B spectral range and most particularly in the UV 
spectral region. In fact, it is mismatch error which 
provides the low temperature cutoff to the spec- 
tral range of stars from which useful UV extinc- 
tion measurements can be made. Mismatch also 
severely limits extinction studies based on stars 
hotter and/or more luminous than the main se- 
quence B stars. For the O stars, very few unred- 
dened standard stars exist and extreme mismatch- 
ing of spectral types is often necessary to derive a 
curve — although, since the intrinsic UV/optical 
SEDs of the O stars are not well known, it is not 
clear how large an effect this introduces in the de- 
rived curves. The use of luminous, evolved stars 
for extinction studies is particularly problematic 
since unreddened standards are rare and the sen- 
sitivity of the intrinsic SEDs to both tempera- 
ture and surface gravity can lead to very severe 
mismatch effects in the resultant curves (although 
see the discussion in Cardelli, Sembach, & Mathis 
1992 for results in the early-B spectral range). 

2.2. Model Atmosphere Techniques 

The second approach to deriving extinction 
curves is to model the intrinsic SED of the red- 
dened star in order to isolate the effects of extinc- 
tion. This technique was first used by Whiteoak 
(1966) to analyze optical spectrophotometry and 
has been applied in one form or another many 
times since. We refer to it as "extinction-without- 
standards", since it does not rely upon a set of 
unreddened standard stars to determine the ex- 
tinction curve. In this method, a model atmo- 
sphere of the reddened star is determined from its 
photometric or spectral properties. The advan- 
tage of this approach is that, in principle, a per- 
fect, unreddened match can be determined for the 
intrinsic SED of the reddened star, eliminating the 
mismatch error which plagues pair method curves. 
The apparent disadvantages of the approach arc 
not actual disadvantages, but rather requirements, 
which can limit its accuracy and range of applica- 
bility. The first requirement is that a set of mod- 
els must exist whose accuracy can be quantified 
and validated. The second is that the observations 



must contain adequate "reddening free" informa- 
tion to accurately determine the intrinsic SED of a 
reddened star. The third is that the fluxes must be 
precisely calibrated, i.e., r\ must well determined. 

We have been developing the necessary con- 
stituents of this method over the past five years. 
We began by demonstrating that the Kurucz 
(1991) ATLAS 9 models provide faithful repre- 
sentations of the observed UV and optical SEDs 
of near main sequence B stars (FM99). We then 
showed that a combination of an observed UV 
SED and optical photometry provides adequate 
reddening independent information to determine 
both the appropriate ATLAS9 model for a red- 
dened star and a set of parameters (defined by 
FM90) which determine the shape of its inter- 
stellar extinction curve. Subsequently, we refined 
the calibration of the IUE data, in order to im- 
prove the quality of the fits and the robustness 
of the physical information derived from them 
(MFOO). Next, we verified the physical parame- 
ters derived from the models through applications 
of the models to eclipsing binary data, where the 
results must agree with other constraints (see, 
Fitzpatrick et al. 2003 and references therein). 
Finally, we have used Hipparcos data of unred- 
dened B stars to derive a consistent recalibration 
of optical and NIR photometry (FM05) and, in 
the process, once again demonstrated the internal 
consistency of the models for near main sequence 
B stars. As a result of these efforts, we now have 
internally consistent r\ for IUE and optical and 
NIR photometry, and are in a position to apply 
the results and to quantify the associated errors. 
These are the objectives of the current paper, and 
are discussed fully in the next section. 

2.3. Normalizing The Curves 

Once color excesses have been determined, we 
are faced with the problem of how to compare 
excesses derived for lines of sight with different 
amounts of extinction. After all, we are inter- 
ested in the "shapes" of the curves, since these 
may reveal important clues about the size distri- 
bution and composition of the dust. This desire 
to compare shapes, brings us to the normalization 
problem. 

From a purely mathematical point of view, a 
straightforward approach would be to normalize 
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the curves by their norm, 

e(X-V)= (9) 

With appropriate weighting, this normalization 
could, on average, minimize the observational er- 
ror in the normalization factor and reduce sys- 
tematic effects. However, such a normalization 
would be sensitive to the strength of narrow fea- 
tures, such as the 2175A bump, and this could 
mask the overall agreement between curves over 
the majority of the wavelength range. A second 
approach is to search for some immutable feature 
of the curve and normalize by that. The idea be- 
hind this procedure is that if some aspect of all 
curves is always the same, then all curves can be 
normalized by the strength of this property and 
all of the resulting curves will be directly compa- 
rable. This is the motivation for the Ay normal- 
ization adopted by Cardelli, Clayton, & Mathis 
(1989). They assume that all extinction curves 
have a very similar (although not identical) form 
for A > V. However, there are problems with this 
normalization as well. In particular, as pointed 
out above, Ay is not a directly measured quan- 
tity, but must be derived from IR photometry and 
requires assumptions about the shape of extinc- 
tion curves at very long wavelengths. The shape 
of this extinction is often considered to have a uni- 
versal form, but this has been demonstrated for 
only a relatively small number of sightlines (Reikc 
& Lebofsky 1985; Martin & Whittet 1990). In ad- 
dition, measurements of extinction in the IR can 
be compromised because the stellar SEDs become 
increasingly faint at long wavelengths and other 
sources of light, such as circumstellar emission or 
scattering from dust in the near stellar environ- 
ment, can contaminate the SED measurements. 
Furthermore, IR color excesses are usually small 
for stars that are detectable in the UV, so UV 
curves normalized by quantities derived from IR 
photometry may be affected by large normaliza- 
tion errors. Consequently, the absolute level of 
such curves can be poorly defined. Finally, even 
with the advent of the 2MASS data base, there are 
still many stars which do not have IR photometry 
available. 

As a result of the complications mentioned 
above, and since our intent is to demonstrate how 
precisely curves can be measured while making a 



minimum number of assumptions, we have opted 
to use the conventional E(B — V) normalization, 
as shown in Eq. 8. While the interpretation of 
E(B — V) as a measure of the "amount" of inter- 
stellar dust is not straightforward, its widespread 
availability, observational precision, and lack of 
requisite assumptions make it the best choice for 
this study. Nevertheless, we note that it is a sim- 
ple matter to transform from one normalization to 
another, and emphasize that E(\ — V) is actually 
the basic measurable quantity. 

3. EXTINCTION- WITHOUT-STANDARDS 

3.1. Formulation of the Problem 

Our earlier studies (FM99, FM05) have shown 
that the observed SEDs, f\, of lightly- or unred- 
dened main sequence B stars can be modeled very 
successfully by using a modified form of Eq. 1, 
namely, 

f\ = F X 6 2 R l()-0-4:E(B-V)[k(\-V)+R(V)] _ ^q) 

The use of absolutely calibrated datasets elimi- 
nates the calibration term r\ and the total ex- 
tinction A\ has been broken down into a normal- 
ized shape term (fc(A — V)), a normalized zero- 
point (R(V) = A V /E(B — V)), and a scale factor 
(E(B — V)). Providing that the righthand side of 
the equation can be represented in a parameter- 
ized form, the equation can be treated as a non- 
linear least squares problem and the optimal val- 
ues of the parameters — which provide the best 
fit to the observed SED fx — can be derived 
along with error estimates. Because the stars un- 
der study were lightly reddened, we could replace 
the extinction curve k(X — V) and the offset term 
R{V) with average Galactic values without loss of 
accuracy. Using the Kurucz ATLAS9 stellar at- 
mosphere models to represent F\, the results of 
the fitting procedure were estimates of 6 parame- 
ters: E(B — V), dji, and the four parameters that 
define the best-fitting model, i.e., T c ff, \ogg, the 
metallicity [m/H], and the microturbulence veloc- 
ity Vt- We performed the fits using the MPFIT 
procedure developed by Craig Markwardt 1 

Markwardt IDL Library available at 

http://astrog.physics.wisc.edu/~craigm/idl/idl.html. 
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The fitting process described above begins to 
break down when the color excess E(B — V) of the 
target stars exceed ~0.05 mag. By this we mean 
that large systematic residuals begin to appear, 
which greatly exceed the measurement errors. The 
reason is simple: the wavelength dependence of 
interstellar extinction curves varies greatly from 
sightline-to-sightlinc, and once E(B — V) > 0.05 
mag the differences between the true shapes of the 
curves and the assumed mean form begin to exceed 
the measurement error. However, FM99 noted 
that the SEDs of significantly reddened stars could 
still be modeled using Equation (10) (and the non- 
linear least squares approach) if the wavelength 
dependence of the extinction curve could be rep- 
resented in a flexible form whose shape could be 
adjusted parametrically to achieve a best fit to the 
observations, and if these parameters were deter- 
mined simultaneously with the stellar parameters. 
This is the essence of our "extinction-without- 
standards" approach. 

Successfully modeling the shape of reddened 
stellar SEDs requires four principal ingredients: 1) 
an observed SED that spans as large a wavelength 
range as possible, 2) an accurate absolute flux cal- 
ibration (r\ or C\), 3) an extinction curve whose 
shape can be described by a manageable set of 
parameters, and 4) a grid of stellar surface fluxes, 
F\, whose defining parameters can be determined 
from the observational data. In §4 we will de- 
scribe the particular datasets used in this paper 
to demonstrate our technique. We have already 
discussed how MF00 and FM05 have determined 
the necessary calibrations. In the remainder of 
this section, we describe our flexible form for the 
interstellar extinction curve (§3.2) and the grid of 
stellar surface fluxes with which we have developed 
our approach (§3.3). 

3.2. A Flexible Representation of the In- 
terstellar Extinction Curve 

We adopt a flexible and adjustable form for the 
UV-through-IR extinction curve, whose shape can 
be optimized to fit the SED of a reddened star 
through the adjustment of a specific set of pa- 
rameters in the least-squares minimization proce- 
dure. This curve is illustrated in Figure 2. It 
consists of two main regions: 1) the UV (A < 2700 
A; solid curve) where the parameterized form of 
FM90 is adopted and 2) the near-UV/optical/IR 



(A > 2700 A; dashed line) where we use a cubic 
spline interpolation through a set of UV (Ui, U2), 
optical (Oi, 2 , 3 , 4 ), and IR (Ii, I 2 , I 3 , I4) 
anchor points to represent the curve. The inter- 
polation is performed using the Interactive Data 
Language (IDL) procedure SPLINE. We adopt a 
spline representation for the near-UV/optical/IR 
curve simply because we do not have reliable, de- 
tailed information on the wavelength dependence 
of the extinction in the near-IR through near-UV 
region (1/im - 3000 A). It is ironic that the por- 
tion of the curve that is accessible from the ground 
is more poorly characterized than the portion ac- 
cessible only from space. As a result, we do not 
know whether the optical to near-IR region of the 
curve can be represented by a compact analytical 
formula. Our hope is that, by applying our proce- 
dure to a large sample of sightlincs, we will ulti- 
mately be able to characterize the shape of the ex- 
tinction law in this region by simple relations and 
determine whether sightline-to-sightline variations 
are correlated with other aspects of the curve or 
with interstellar environment. The placement of 
the spline points resulted from considerable exper- 
imentation, but certainly cannot be represented 
as an objectively determined optimal result. The 
current arrangement does, however, allow us to 
model the major available datasets to a level con- 
sistent with the observational errors. 

The FM90 parameterization scheme contains 6 
free parameters to represent 3 functionally sepa- 
rate features that are summed to produce the UV 
curve. An underlying linear component, indicated 
by the dotted line in Figure 2, is specified by an 
intercept c\ and a slope C2. The Lorentzian-likc 
2175 A bump is fit by a Drude profile D(x,x ,j), 
where x$ and 7 specify the position and FWHM 
of the bump, respectively, whose strength is deter- 
mined by a scale factor C3. Finally, the degree of 
departure of the curve in the far-UV from the un- 
derlying linear component is specified by a single 
parameter C4. Defining x = A -1 , the complete UV 
function is given by: 

fc(A - V) = ci + c 2 x + c 3 D(x, x , 7) + c 4 F(x) , 

(11) 

where 

^^Tl^^^, (12) 
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and 

F(x;x> 5.9) = 0.5392(a; - 5.9) 2 + 0.05644(x - 5.9j(l3) 
F(x;x < 5.9) = . (14) 

It is worth emphasizing that the FM90 param- 
eterization is a mathematical scheme only, which 
allows us to reproduce UV extinction curves in a 
shorthand form (and makes the current program 
possible). However, it is not to be assumed that 
the functional components of the scheme represent 
actual separate extinction components arising in 
distinct dust grain populations. This parameter- 
ization has proven very useful and to the best of 
our knowledge is able to reproduce all known UV 
extinction curves to the level of observational er- 
ror. Its flexibility will be illustrated in §5 below. 
Note that we terminate the FM90 formula at 2700 
A, although we originally utilized UV data extend- 
ing to 3000 A. Additional experience with the UV 
data has suggested that real extinction curves be- 
gin to exhibit departures from the FM90 formula 
near 2700 A, due to the unrealistic extrapolation 
of the linear component into the blue- violet region. 

The 10 spline anchor points which character- 
ize the near-UV/optical/IR portion of the curve 
are determined by a least squares fit to the IUE 
data longward of 2700 A and the available optical 
and IR photometry. Although there are 10 an- 
chor points, there are also 6 constraints, so the 
fit actually introduces only 4 additional degrees 
of freedom. We now discuss these constraints in 
detail. 

The two UV anchor points, Ui and U2 at 2700 
and 2600 A, respectively, are fixed at the values 
of the FM90 UV fitting function at their wave- 
lengths and are not adjustable. Together with Oi 
at 3300A, these points guarantee that the curve 
which passes through the IUE data between 2700 
and 3000 A will join both the UV and optical por- 
tions of the curve smoothly. 

The four optical anchor points, Oi, O2, O3, and 
4 at 3300, 4000, 5530, and 7000 A, respectively, 
are fit under two constraints: that the interpolated 
curve produces a value of k(X — V) = in the V 
band, and that the curve be normalized to unity in 
E(B — V). Thus, only two free parameters emerge 
from this region. 

The four IR points, T, I2, I3, and I4 are located 
at 0.25, 0.50, 0.75, and 1.0 fim -1 , respectively. 
These four points are constrained to satisfy the 



formula 

I n = k(\ -V) = k IR \- 1M - R(V) (15) 

where the scale factor, km, and the intercept, 
R(V), are the only free parameters. This is the 
power-law form usually attributed to IR extinc- 
tion, with a value for its exponent from Martin & 
Whittet (1990). The exponent of the power-law 
could, potentially, be included as a free param- 
eter in the fitting procedure, and we will investi- 
gate this in the future. However, our impression is 
that the IR data available to us (primarily 2MASS 
JHK photometry) are insufficient to determine 
this quantity accurately. All results presented in 
this paper assume an exponent of -1.84 in Equa- 
tion 15. 

3.3. The Stellar Surface Fluxes 

To represent the intrinsic surface fluxes, F\, 
of reddened stars wc utilize R.L. Kurucz's line- 
blanketed, hydrostatic, LTE, plane-parallel AT- 
LAS9 models, computed in units of erg cm -2 sec" 1 A 
and the synthetic photometry derived from the 
models by Fitzpatrick & Massa (2005). These 
models are functions of four parameters: T e g, 
logy, [m/H], and vt- All of these parameters can 
be determined in the fitting process although, be- 
cause of data quality, it is sometimes necessary to 
constrain one or more to a reasonable value and 
solve for the others. 

The general technique of deriving extinction 
curves via stellar atmospheres is, of course, not 
dependent on the specific set of models used. In 
the present case, the most important considera- 
tion for our adoption of the ATLAS9 models is 
that they work — at least within a specific spec- 
tral domain. FM99 and FM05 have shown that 
these models provide excellent fits to the observed 
SEDs for lightly- or un-reddened main sequence 
B stars throughout the UV through near-IR spec- 
tral region. In addition, experience with eclipsing 
binaries (see Fitzpatrick et al. 2003 and refer- 
ences within) has shown that the good SED fits 
also yield accurate estimates of the physical prop- 
erties of the stars. Because of the physical ingre- 
dients of the models (specifically LTE and plane- 
parallelism) we currently restrict our attention to 
the main sequence B stars. We plan to investigate 
how well these models reproduce the SEDs and 
properties of somewhat more luminous B stars and 
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also the later O types. Also, we will take advan- 
tage of more complex models (e.g., the non-LTE 
TLUSTY models) as the available grids expand 
their parameter ranges. 

3.4. Summary 

To summarize, we model the observed SEDs of 
reddened near main sequence B stars by treating 
equation (10) as a non-linear least squares mini- 
mization problem. As a result we can simultane- 
ously obtain estimates of the physical properties of 
a reddened star and the shape of the interstellar 
extinction curve distorting the star's SED. A total 
of 16 parameters specify the righthand side of the 
equation, including Or, E(B — V), four to define 
F\, and ten to define the shape of the extinction 
curve. Depending on limitations of the available 
data, and known properties of the stars or inter- 
stellar medium, any subset of the parameters can 
be constrained to predetermined values. 

4. THE DATA 

In the following section, we will illustrate 
the potential of our extinction-without-standards 
technique, utilizing a set of 27 lightly-to-heavily 
reddened stars. For this demonstration, and in- 
deed for extinction determinations in general, the 
ideal SED datasct would consist of absolutely- 
calibrated spectrophotometry spanning the UV- 
through-IR spectral regions. Such data would al- 
low a straightforward comparison between obser- 
vations and stellar atmosphere models (since both 
are presented in simple flux units) and would pro- 
vide the most detailed view of the wavelength de- 
pendence of interstellar extinction. While a small 
amount of such data is available (see, e.g., Fitz- 
patrick at al. 2003), the largest existing database 
of absolutely-calibrated spectrophotometry is the 
low-resolution archive of the International Ultra- 
violet Explorer satellite (IUE) which covers the 
UV region only (1150-3000 A). In the optical and 
near-IR, the largest SED databases are photomet- 
ric in nature, consisting of Johnson, Stromgrcn, 
and Geneva photometry in the optical and 2MASS 
JHK photometry in the near IR. Utilizing these 
resources, we can examine the UV region for small 
scale features but can only study the broad scale 
wavelength-dependence of extinction in the opti- 
cal and near-IR regions. 



We use NEWSIPS IUE data (Nichols & Linsky 
1996) obtained from the MAST archive at STScI. 
These data were corrected for residual systematic 
errors and placed onto the HST/FOS flux scale 
of Bohlin (1996) using the corrections and algo- 
rithms described by MF00. This step is abso- 
lutely essential for our program since our "com- 
parison stars" are stellar atmosphere models and 
systematic errors in the absolute calibration of the 
data do not cancel out as in the case of the pair 
method. (The NEWSIPS database is also contam- 
inated by thermally- and temporally-dependent 
errors, which would not generally cancel out in the 
pair method — see MF00.) Multiple spectra from 
each of lUE's wavelength ranges (SWP or LWR 
and LWP) were combined using the NEWSIPS er- 
ror arrays as weights. Small aperture data were 
scaled to the large aperture data and both trailed 
and point source data were included. Short and 
long wavelength data were joined at 1978 A to 
form a complete spectrum covering the wavelength 
range 1150 < A < 3000 A. Data longward of 
3000 A were ignored because they are typically 
of low quality and subject to residual systematic 
effects. The IUE data were resampled to match 
the wavelength binning of the ATLAS9 model at- 
mosphere calculations in the wavelength regions of 
interest. 

Mean values of the Johnson UBV, Stromgrcn 
uvbyfi, and Geneva UBB\BiVV\G photometric 
magnitudes, colors, and indices for the program 
stars were acquired from the Mermilliod et al. 
(1997) archive. 2MASS JHK magnitudes for all 
stars, along with their associated errors, were ob- 
tained from the 2MASS All-Sky Point Source Cat- 
alog at the NASA/IPAC Infrared Science Archive. 
Johnson V—R, R—I, and JHK data, are also avail- 
able for some of the stars, and were obtained from 
the Mermilliod et al. archive. 

5. SOME INITIAL RESULTS 

In this section, we demonstrate the potential of 
our extinction-without-standards technique, uti- 
lizing a set of the 27 reddened stars, listed in Ta- 
ble 1, which fall into three groups: 1) stars in the 
open cluster IC 4665; 2) stars with moderate-to- 
heavy reddening; and 3) lightly reddened stars in 
a specific region of the sky. These representative 
examples illustrate the advantages of our approach 
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and highlight several scientific applications which 
will be pursed in future studies, using expanded 
samples of stars. In addition, they provide con- 
firmation of the error analysis incorporated in our 
approach. 

For this demonstration sample, the SED fitting 
procedure was applied as described above, with 
the following additional details: 

• The SED data modeled in the fitting pro- 
cedure include IUE UV spectrophotometry, 
the Johnson V, B — V, and U — B indices, 
the Stromgren b — y, mi, Ci, and indices, 
the Geneva U - B, V - B, B x - B, B 2 - B, 
Vi - B, and G - B indices, and the 2MASS 
JHK magnitudes. Johnson V — R, R — I, 
and JHK data are also available for a few of 
the stars. 

• The optical extinction spline point O4 at 
7000 A is only well-determined when opti- 
cal R and / band photometry are available. 
Therefore in the examples below, we only 
solve for O4 in such cases. For the other 
stars, the optical portion of the extinction 
curve is determined only by the spline points 
Oi, 2 , and 3 . 

• During the x 2 minimization, a reddened and 
distance-attenuated model was created from 
the current set of input parameters, and 
then synthetic photometry was performed 
on this model to produce the photometric 
indices, which were then compared with ob- 
servations. The synthetic photometry was 
calibrated as described by FM05, with the 
calibration extended to redder colors by us. 
The UV model fluxes and recalibrated IUE 
fluxes, both in units of erg cm -2 sec -1 A -1 , 
were compared directly. 

• The initial weighting factors for the vari- 
ous datasets in the \ 2 minimization were 
determined from their observational uncer- 
tainties, i.e., weight oc a~ 2 . We then scaled 
the weights of the optical/near-IR photom- 
etry so that their total weight was equal to 
that of the IUE UV spectrophotometry, thus 
balancing the fit between the two datasets. 
This is the procedure adopted by FM05, ex- 
cept that we include the Stromgren index 
along with the rest of the optical/near-IR 



photometry, rather than assigning it its own 
(high) weight. FM05 weighted heavily in 
recognition of its value as a surface gravity 
indicator. However, we have found that the 
temperature-sensitivity of (particularly in 
the later B stars) combined with observa- 
tional errors, can lead to very unsatisfactory 
fits when is over-emphasized. Treating 
in the same manner as the rest of the photo- 
metric indices seems to be the simplest and 
most reasonable approach. 

• Because interstellar H 1 Lyman a absorp- 
tion in reddened stars can have a signifi- 
cant impact on the star's apparent contin- 
uum level far from line center at 1215 A, we 
convolve the profile of this heavily damped 
line with the model atmosphere SEDs be- 
fore comparing with observations. Along 
heavily reddened sightlines, where the H 1 
column density N(H 1) is high and the sig- 
nature of the atomic absorption strong, the 
value of N(R 1) can actually be incorporated 
into the fitting procedure as a free parame- 
ter to be optimized. We will utilize this ca- 
pability in future studies. For the present, 
mostly lightly-reddened, sample, however, 
the iV(H 1) values used to construct the 
line profiles were taken from the survey of 
Diplas & Savage (1994) or else computed 
from the general relation N(H 1) = 4.8 x 
10 21 E(B - V) cm" 2 from Bohlin, Savage, & 
Drake (1978). The inclusion of the Lyman 
a line insures that we distinguish the effects 
of dust extinction from atomic absorption in 
the far UV region. 

• The uncertainties in the best-fit parame- 
ters were determined by running 50 Monte 
Carlo simulations for each star, during which 
the input data were randomly varied as- 
suming a Gaussian distribution of obser- 
vational uncertainties and a new fit per- 
formed. The zero points and random photo- 
metric uncertainties of the short-wavelength 
and long-wavelength IUE fluxes were varied 
as described in FM04; the assumed obser- 
vational errors in the Johnson, Stromgren, 
and Geneva indices were as given in Ta- 
ble 7 of FM04; and the uncertainties in 
the 2MASS data were as obtained from the 
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2MASS archive. In addition, the V magni- 
tude was assumed to have a l-er uncertainty 
of 0.015 mag. The adopted 1-a uncertainties 
for each parameter were taken as the stan- 
dard deviation of the values produced by the 
50 simulation. 

5.1. The Open Cluster IC 4665 

We begin our examples by examining the ex- 
tinction towards an open cluster. While multiple 
scientific rewards can result from the study of ex- 
tinction towards cluster stars (see the discussion in 
§6), our primary interest in IC 4665 is to demon- 
strate the "technical" advantages of our approach. 
Namely, the use of cluster extinction curves to help 
evaluate the magnitude of the uncertainties in the 
measurement of extinction curves, as discussed in 
detail by Massa & Fitzpatrick (1986). In particu- 
lar, and because of its low E(B — V) and late-B 
stellar population, IC 4665 extinction curves pro- 
vide an especially sensitive test of the precision 
and range of our extinction- without-standards ap- 
proach. 

The wavelength dependence of extinction to- 
wards IC 4655 was first examined by Hackwell, 
Hecht, & Tapia (1991; hereafter HHT) for the pur- 
pose of studying the relationship between extinc- 
tion and IR emission, as measured by the Infrared 
Astronomical Satellite (IRAS). This remains one 
of the most challenging extinction studies yet per- 
formed for two reasons: 1) the mean reddening in 
the cluster is very low, (E(B — V)) < 0.2 mag, 
and 2) the spectral types of the available targets 
run from mid- to late B. Both facts exacerbate er- 
rors in the standard pair method approach, as has 
been shown in Figure 1. HHT recognized these 
uncertainties and ultimately concluded that the 
wavelength dependence of extinction among the 
cluster stars is uniform to within their ability to 
measure it. 

Figure 3 shows the results of the SED fits for the 
nine IC 4665 stars considered here. The SED's of 
the best-fitting, reddened models are shown by the 
histogram-style curves. In the UV, the binned IUE 
fluxes are shown by the small circles. In the optical 
region, Johnson UBVRI magnitudes (converted to 
flux) are indicated by circles, Stromgren uvby mag- 
nitudes by triangles, and Geneva UBiB 2 VViG 
magnitudes by diamonds. In the near-IR, 2MASS 
and Johnson JHK magnitudes are shown by the 



large filled and open circles, respectively. Note 
that the photometric data have been converted to 
flux units for display purposes only. The com- 
parison between models and observations was per- 
formed in the native photometric format (i.e., in 
magnitudes or colors as noted above). 

Figure 4 shows our extinction-without-standards 
curves for the IC 4665 stars. The symbols show 
the actual normalized ratios between the models 
and the stellar SEDs, while the thick solid curves 
show the flexible UV-through-IR extinction curves 
whose shapes were determined by the fitting pro- 
cedure. The curves have been arbitrarily shifted 
vertically for clarity, but all are shown compared 
with a similarly-shifted estimate of the average 
Galactic extinction curve for reference (thick dash- 
dot curves corresponding to R{V) = 3.1, from 
Fitzpatrick 1999). As will be discussed further 
below, we assumed a value of R(V) =3.1 towards 
the cluster and did not include the IR JHK data 
in the fitting procedure. Thus, only the average 
Galactic curve is shown for wavelengths longward 
of 6000 A. 

The various parameters determined by the fits 
are listed in Tables 1 and 2. The flexible extinction 
curves themselves, in the form E(\ — V)/ E(B — V) 
can be reconstructed from the parameters given 
in Table 2. The l-c uncertainties of the extinc- 
tion curves are indicated in Figure 4 by the grey 
shaded regions. The regions are centered on the 
means of the 50 Monte Carlo simulations with 
which we performed our error analysis and their 
thickness shows the standard deviation of the in- 
dividual simulations. 

Figure 5 compares our new curves in the UV 
with those derived by HHT using the standard 
pair method. HHT's curves were reconstructed 
from the data in their Table 5. Note that the 
HHT study originally included 17 stars. We have 
eliminated 4 A-type stars, 3 chemically peculiar B- 
type stars, and 1 B-type shell system from consid- 
eration here since their extinction curves are par- 
ticularly uncertain. Thus, Figures 4 and 5 show 
only the best-determined curves in the cluster, 
from the point of view of both the pair method 
and our technique. The remarkable aspect of 
Figure 5 is not the scatter among HHT's curves 
- it is exactly what should be expected given 
the limitations of the pair method — but rather 
the impressive improvement in precision gained by 
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the extinction- without-standards technique, as ev- 
idenced by the decrease in the curve-to-curve scat- 
ter. Clearly, this higher level of precision allows 
much firmer conclusions to be drawn about the 
degree of intrinsic variation of extinction across 
the face of the cluster, and also affords an im- 
proved potential to detect correlated behavior be- 
tween the extinction properties and other aspects 
of the interstellar environment. 

While we will present a full scientific analysis 
of the IC 4665 results in a future paper, the issue 
of the intrinsic variation of extinction among the 
cluster stars is important to consider here, since it 
can shed light on the accuracy of our error analy- 
sis. The top curve in Figure 4 shows the mean of 
the nine individual IC 4665 curves. The error bars 
show the actual sample standard deviation and the 
average Galactic curve is again presented for com- 
parison. The critical result is that, over most of 
its wavelength range, the standard deviation of the 
sample about the mean curve is comparable to the 
predicted uncertainties in the individual curves, as 
shown by the shaded regions. This indicates that 
the (small) curve-to-curve variations seen are at 
a level consistent with our expected uncertainties 
and that the intrinsic level of variation among the 
cluster stars must be very low. 

Another way to approach the issue of variability 
is to look at the scatter among the various param- 
eters which define the extinction curves. These 
are shown in Table 3, where columns 2 and 3 list 
the weighted mean values of the parameters and 
their observed standard deviations, respectively. 
The predicted uncertainties, i.e., the RMS of the 
Monto Carlo-based errors for the individual stars, 
arc listed in column 4. If our error analysis is 
reasonable, then the observed scatter should be 
the quadratic sum of the expected errors plus any 
intrinsic variability. The value of examining clus- 
ter extinction curves is that one might reasonably 
suppose (at least as a starting point) that the in- 
dividual curves, derived for nearly coincident lines 
of sight, are actually independent measurements 
of a single "cluster curve," i.e., no intrinsic scat- 
ter. In such a case, the measured scatter actually 
reflects the real measurement errors and provides 
an important test of the error analysis. Examina- 
tion of Table 3 shows that the predicted and ob- 
served scatters are indeed generally very similar, 
supporting the position that any intrinsic varia- 



tions among the IC 4665 curves are close to the 
level of our ability to measure them and that we 
have accurately assessed the uncertainties in our 
results. The final column of Table 3 shows the 
implied values of the possible small intrinsic vari- 
ations. 

Several of the individual parameters merit some 
additional comment. Both the 2175 A bump 
FWHM (7) and its strength (c 3 or Ab ump ) show 
evidence for some small level of variability within 
the sample, above our expected measurement er- 
rors. However, this is not clearcut because these 
measurements involve the region of the IUE spec- 
tra which typically has the lowest quality data and 
it is possible that weak systematic effects in the 
data themselves — which are not accounted for 
in the error analysis — could produce the small 
level of variability seen. One such systematic is a 
"reciprocity effect" in long wavelength data which 
is not corrected by the MFOO algorithms (see Fig- 
ures 12 and 13 of MFOO). We conclude conser- 
vatively that there is marginal evidence for bump 
variations within the IC 4665 sample, but will ul- 
timately rely on studies of several open clusters 
to determine whether our error analysis faithfully 
predict the real measurement errors in the bump 
region. 

The case of the far-UV curvature (parameter 
C4), for which the implied intrinsic scatter is three 
times greater than our measurement errors, is 
more interesting. It is clear from Figures 4 and 5 
and Table 3 that most of this apparent variability 
arises from two sightlines, namely, those towards 
HD 161165 and HD 161184. In fact, the observed 
scatter in C4 towards the other seven sightlines is 
essentially identical to the expected measurement 
error. Because HD 161165 and HD 161184 arc the 
two coolest stars in the sample, and the only stars 
with effective temperatures less than 12,000 K, we 
suspected that the high C4 values in their extinc- 
tion curves might be artifacts of the analysis, re- 
sulting from a failure of the models to accurately 
portray the far-UV SEDs of these stars. The inves- 
tigation presented below suggests that this could 
well be the case, with the extinction curve anoma- 
lies possibly arising from a difference between the 
chemical composition profiles of the stars and the 
atmosphere models. 

The data in Table 1 show a very uniform 
"metallicity" [m/H] for the cluster, with a weighted 
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mean of -0.50 and a standard deviation of 0.04. 
This is close to, and actually slightly smaller than, 
our estimate of the measurement errors, simulta- 
neously confirming the accuracy of our error anal- 
ysis and imposing a small upper limit on the in- 
trinsic compositional variability within the cluster. 
Although our [m/H] values are simple scale factors 
which apply to a template set of ATLAS9 solar 
abundances, FM99 showed that the [m/H] derived 
in our analysis is most sensitive to the abundance 
of Fe — due to a very strong opacity signature in 
the mid-UV — and is most analogous to [Fe/H]. 
Since the Fe abundance in the ATLAS9 models 
is ~0.2 dex larger than the currently accepted 
solar value of 7.45 (where H = 12.00; Asplund, 
Grevesse, & Sauval 2005), our results suggest that 
the IC 4665 stars are deficient by about a factor 
of two in Fe as compared to the Sun. The small 
scatter observed in [m/H] is both satisfying and 
expected, since the surface composition of Fe is 
not subject to evolutionary modification in these 
young stars, which presumably all formed from 
the same parent material. 

We experimented with forcing [m/H] to a more 
solar-like level for the IC 4665 stars and imme- 
diately found two effects: 1) the x 2 values all 
increased significantly because the Fe features in 
the mid-UV could not be fit as well, and 2) most 
of the extinction curves remained unchanged, but 
the anomaly in the far-UV region for HD 161165 
and HD 161184 decreased dramatically. The first 
result was expected. The second was a surprise, 
but is understandable. If the cluster stars (or at 
minimum, the two coolest stars HD 161165 and 
HD 161184) have a non-solar ratio of Fe to the 
light metals, e.g., [C/Fe] > 0, then our best fit 
models — biased towards the Fe abundance - 
would underestimate the opacity due to the light 
metals. For the cooler stars, such elements provide 
significant opacity in the far-UV region and the 
ATLAS9 models would not be able to account for 
such a specific opacity difference. The fitting pro- 
cedure would respond by finding a higher far-UV 
extinction curve. The curves for the hotter stars 
are less affected by changing [m/H] since the light 
metal opacity is less significant. We tentatively 
conclude that the chemical composition of the B 
stars in IC 4665 may deviate strongly from that 
of the Sun, with subsolar Fe but a more "normal" 
level of the light metals such as C. This suggestion 



is easily tested with a fine analysis of high resolu- 
tion stellar spectra and we will pursue this in the 
future. On the positive side, this result suggests 
that the UV continua of late B and early A stars 
might be used to determine both a scaled [m/H] 
and a mean [CNO/Fe] index, assuming a grid of 
models parameterized by both these composition 
indices is available. On the more sober side, it is a 
reminder that our technique is susceptible to stel- 
lar abnormalities. As with the pair method, all 
available data should be consulted to determine 
whether a particular star is suitable for deriving 
an extinction curve. 

Our final comment on the IC 4665 results con- 
cerns the behavior of the IR data. In Figure 4 we 
show two sets of JHK photometry for each star 
and for the average cluster curve. The solid sym- 
bols indicate 2MASS measurements and the open 
symbols show HHT's data. The two sets of data 
systematically differ, with the 2MASS results sug- 
gesting a value of R(V) somewhat smaller then the 
average Galactic value of 3.1 and the HHT data 
suggesting a slightly larger value (HHT derive a 
mean of R{V) = 3.25). This latter result is more 
consistent with expectations, given that the mean 
far-UV curve is lower than the Galactic average 
(see, e.g., Cardelli et al. 1989), but this is in- 
sufficient evidence for rejecting the 2MASS data. 
For now, our solution for this quandary has been 
to ignore both sets of data in fitting the IC 4665 
SEDs and adopt a default value of R(V) = 3.1. 
However, the discrepancy for measurements in this 
very complex region bears further investigation, 
as does that fact that, in both datasets, the mean 
curves actually show more extinction in the 2.2 /xm 
K band than in the 1.65 /zm H band. 

The discussion above leads us to three primary 
conclusions: 1) The extinction towards IC 4665 
shows at most only a small degree of spatial vari- 
ability comparable with our measurement er- 
rors — and that the cluster extinction curve would 
best be represented by averaging the results for 
the seven hottest cluster stars studied here; 2) 
we are able to determine accurate [m/H] values 
from the observations; and 3) our analysis yields 
reliable estimates of the (small) uncertainties in 
the extinction-without-standards curves and pa- 
rameters, barring the presence of unusual system- 
atic anomalies in the data sample. The first and 
second conclusions are scientific issues which we 
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will pursue further in the future. The third pro- 
vides a formidable demonstration of the superior- 
ity of the extinction-without-standards technique 
over the classical pair method for deriving extinc- 
tion curves in both the precision of the results and 
the quantification of the uncertainties. 

5.2. Moderately- Reddened Stars 

With our error estimates verified, we now exam- 
ine how curves derived by the current approach 
compare to pair method curves. Figure 6 shows 
the UV-through-IR fits to the SEDs of a set of 9 
moderately-reddened early-B stars, most of which 
are well known for their extinction properties, 
and Figure 7 shows the corresponding extinction 
curves. These stars were selected because they 
illustrate the wide range that exists in the wave- 
length dependence of both UV and IR extinction. 
For these stars the IR data allow us to determine 
the values of R(V) and so the flexible extinction 
curve fits are shown throughout the IR-to-UV do- 
main. As for the IC 4665 stars, all the parameters 
describing the fits are given in Tables 1 and 2. 

Also shown in Figure 7 are UV curves based 
on the pair method technique. The curve for 
HD 294264 is from Valcncic, Clayton, & Gordon 
(2004); the curves for HD 210121 and HD 27778 
are from unpublished measurements by us; and the 
others are from the catalog of FM90. The agree- 
ment between the model-based and pair method 
curves is reasonable, and much better than seen for 
IC 4665. This is consistent with spectral mismatch 
as the prime cause of the existing discrepancies, 
since the nine stars have both higher E(B — V) 
values and earlier spectral types than the IC 4665 
stars, both of which tend to reduce the influence 
of mismatch errors. Note also that the best agree- 
ment between the pair method and the model- 
based curves occurs for the star HD 204827, for 
which identical results are found. Again, this is as 
it should be, since its E(B — V) is the largest of 
any star in the sample, minimizing the impact mis- 
match effects. The 1-a uncertainties for our flex- 
ible extinction curve fits are indicated by shaded 
regions around the curves, as in Figure 4. 

The curves shown in Figure 7 demonstrate the 
ability of the parameterized, flexible extinction 
curve (see §3.2) to conform itself to the wide range 
of extinction curve shapes encountered in interstel- 
lar space. 



5.3. Lightly-Reddened Stars 

Figure 8 shows a set of nine SED fits for low 
color excess stars (0.10 < E(B — V) < 0.21) lo- 
cated along sightlines bounded by the Galactic co- 
ordinates 347° < I < 355° and 18° < b < 26°. Fig- 
ure 9 shows the corresponding extinction curves. 
These stars are all mid-to-late B members of the 
Upper Scorpius complex (Garrison 1967) and we 
encountered them while testing a procedure for 
scanning the IUE archives and automatically gen- 
crating model-based extinction curves for stars in 
the appropriate spectral range. When examining 
the results for this pilot program which sampled 
high-latitude stars, it became obvious that curves 
derived from stars in this specific region demon- 
strated similar curve morphology which is dis- 
tinctly different from the Galactic average. Given 
the low reddening and preponderance of late B 
stars, this strong systematic behavior would be 
missed by a pair method survey, with its large in- 
herent mismatch errors. 

Strong regional signatures are important in ex- 
tinction studies because they may highlight the 
effects of specific physical processes on dust grain 
populations. Although we will examine this spe- 
cific region in the future, two points are worth 
mentioning. First, it is important that similar 
curves result from stars with spectral types rang- 
ing from B2.5 to B9, verifying that the curves are 
not the result of some temperature-dependent sys- 
tematic effect in fitting process. It is also worth 
noting that the region is near the p Oph dark 
cloud. The star p Oph A (HD 147933, whose ex- 
tinction curve is shown in Figure 7) is located at a 
comparable distance and just south of the region, 
at Galactic coordinates of (I, b) = (353.7°, 17.7°). 
The mean curve for the nine lightly-reddened stars 
is shown at the top of Figure 9, along with the 
average Galactic curve (dash-dot curve) and the 
HD 1479433 curve (dotted curve) for comparison. 
In the UV, the regional curve is seen to be almost 
identical to that for the more heavily-reddened 
HD 147933 (E(B - V) ~ 0.5), suggesting that 
the sightlines sample similar dust populations. In- 
terestingly, however, the curves are not identical 
in the IR. The mean 2MASS JHK data for the 
nine star sample implies a value of R(V) ~ 3.4, 
slightly larger than the Galactic mean of 3.1, while 
HD 147933 has a value of R(V) ~ 4.3. This 
contrast between UV and IR behavior may in- 
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dicate different physical processes at work along 
the higher density sightline towards HD 147933, or 
perhaps different timescales in the response of UV 
and IR extinction to modifications of dust grain 
properties. 

In producing the extinction curves for these 
lightly reddened stars, we found that, in some 
cases, the value of km was very poorly determined 
and produced (presumably) spurious "bumps and 
wiggles" in the IR portion of the curves. To elimi- 
nate this effect, we imposed a constraint on km for 
the whole sample, namely km = 0.63R(V) — 0.84. 
This is taken from F04, who found a very strong 
relationship between R(V) and km from a larger 
sample of more heavily reddened stars (see Fig- 
ure 6 of F04). The ability to impose scientif- 
ically reasonable constraints on the fitting pro- 
cedure is a major advantage of the extinction- 
without-standards approach, and can potentially 
allow high quality extinction results to be derived 
from stars with even lower reddenings than those 
shown in Figure 9. 

6. DISCUSSION 

In the previous sections we first provided an 
overview of the process used to determine an ex- 
tinction curve, clarifying the measurement pro- 
cess. We then presented a new method for de- 
riving UV-through-IR extinction curves, based 
on the use of stellar atmosphere models to pro- 
vide estimates of the intrinsic (unreddened) stel- 
lar SEDs rather than unreddened (or lightly red- 
dened) "standard" stars. We have shown that this 
"extinction-without-standards" technique greatly 
increases the accuracy of the derived extinction 
curves, particularly in the cases of low reddening 
and cool spectral types (i.e., late-B), and allows 
a realistic estimation of the uncertainties. A side 
benefit of the technique is the simultaneous de- 
termination of fundamental properties of the red- 
dened stars themselves (T c g, log g, [m/H], and 
Vturb), making the procedure valuable for both 
stellar and interstellar studies. Given the physical 
limitations of the ATLAS9 models we currently 
employ, the technique is limited to near main se- 
quence B stars. However, in principal, the proce- 
dure can be adapted to any class of star for which 
accurate model SEDs are available and for which 
the signature of interstellar reddening can be dis- 



tinguished from those of the stellar parameters. 
Although we developed the procedure based on 
IUE spectrophotometry in the UV, and photome- 
try in the optical and near-IR (requiring a calibra- 
tion of synthetic optical and near-IR photometry), 
the ideal application of the technique would be 
with spcctrophotomctric data throughout the UV- 
through-IR domain, allowing the most detailed ex- 
amination of the wavelength dependence of the ex- 
tinction curves. 

The specific scientific advantages afforded by 
the extinction-without-standards technique can be 
summarized as follows: 

Increased Precision: The increased precision in 
the derived extinction curves allows us to improve 
our understanding of extinction in a number of 
ways. First and most simply, we will determine 
the basic UV-to-IR wavelength dependence of ex- 
tinction along a wide variety of sightlines more 
precisely than has been possible in the past. Sec- 
ond — and given that we know that extinction 
curve morphology varies widely from sightline-to- 
sightline (see Figure 7!) — we will be able to study 
the form of the variability and search for relation- 
ships between various features and wavelength do- 
mains using a data set with small and well-defined 
uncertainties. Such relationships, e.g., the corre- 
lation between R(V) and the steepness of UV ex- 
tinction discovered by Cardelli et al. (1989), pro- 
vide important constraints on the dust grain pop- 
ulation causing the extinction. Non-correlations 
can be equally important. For example, the lack 
of a correlation between the position and width 
of the 2175 A bump demonstrated by Fitzpatrick 
& Massa (1986) remains a strong constraint on 
models for the bump carrier (Draine 2003). In ci- 
ther case, a precise knowledge of the measurement 
errors is required for transforming an observation 
into a scientific constraint. Finally, we will be able 
to place much stronger and more well-defined lim- 
its on the relationship between extinction curve 
morphology and interstellar environment. Curves 
derived from lines of sight to specific, localized re- 
gions, such as towards open clusters, can be partic- 
ularly useful for relating curve properties to phys- 
ical processes occurring in the interstellar medium 
(see, e.g., the study of Cepheus OB3 by Massa & 
Savage 1984 and Trumpler 37 by Clayton & Fitz- 
patrick 1987). 

Access to lightly-reddened sightlines: The 
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ability to accurately probe low E{B — V) sightlines 
(as exemplified in Figures 4 and 9) opens the door 
to studies of dust in regions that have not been 
thoroughly explored. These include halo dust, 
dust in very low density regions, and local dust. 
Halo dust is especially important since we must 
contend with its effects every time wc look out 
of the Galaxy. There are indications (Kiszkurno- 
Koziej & Lequeux 1987) that its properties differ 
systematically from dust at lower Galactic alti- 
tudes, and this result needs to be verified on a 
star-by-star basis. The nature of low density dust 
provides insights into the processing which occurs 
in hostile environments. Clayton et al. (2000) 
presented observations of dust from low density 
sightlines, and their results were intriguing. How- 
ever, they were forced to select sightlines which 
accumulate fairly substantial color excesses, intro- 
ducing the possibility of mixed grain populations. 
Furthermore, the results for their least-reddened 
sightlines were, as they acknowledge, poorly de- 
termined, forcing them to base their conclusions 
on a global average of properties of their sample. 
Finally, measuring the properties of local dust is 
important because it allows us to search for iso- 
lated, relatively homogeneous environments with 
uniform extinction curve shapes. These may sig- 
nal physically and kinematically isolated regions 
which would be ideal for follow up interstellar line 
studies. In addition, Hipparcos parallax data ex- 
ist for nearby stars and provide an opportunity to 
study the 3-dimensional structure of local extinc- 
tion. 

Access to the mid-to-late B stars: These stars 
are especially important because their space den- 
sity is higher than that of the early-B stars usually 
used in extinction studies and thus their inclusion 
increases the number of stars available to cre- 
ate curves for nearby sightlines. Study of the 
mid-to-late B stars will enable the examination 
of the spatial structure of local dust absorption 
more thoroughly than previously possible and will 
greatly enhance our understanding of the local 
interstellar medium. The ability to construct ac- 
curate curves for mid-to-late B stars will also allow 
us to study extinction in open clusters and asso- 
ciations which are dominated by these stars, such 
as the Pleiades, a Per cluster, and IC 4665 (see 
Figure 4). 

Automation: Because our model atmosphere ap- 



proach does not require human intervention, once 
the basic data have been assembled, it is possible 
to process large data sets at one time. Naturally, 
the automated results must be inspected for out- 
liers and data anomalies. Nevertheless, this ap- 
proach reduces the work load considerably and a 
first-attempt yielded the regional anomaly shown 
in Figure 9. Furthermore, since the results are 
produced in a uniform manner, it is a relatively 
simple matter to inspect them for correlations 
between various curve properties, for anomalous 
curve shapes, and for spatial trends on the sky. 
Stellar properties: In addition to dust parame- 
ters, our technique provides a meaningful, quanti- 
tative physical properties for the reddened stars. 
The temperature and surface gravity information 
will be useful for population studies of B stars in 
the field and in clusters. However, perhaps the 
most useful property will be the metallicity. We 
have demonstrated that several stars in the same 
cluster, which have a range in temperatures and 
gravities, all yield the same [m/H]. This verifies 
the sensitivity of our fitting procedure to this im- 
portant quantity. As a result, we are confident 
that application of our procedure to large scale 
surveys of reddened, near main sequence B stars 
can provide a census of the distribution of metal- 
licity throughout the Galaxy and the local uni- 
verse. 
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Table 1 

Best-Fit Parameters for Stars in Figures 3-9 
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a The maximum allowed value of the surface gravity was constrained to be (log g) m ax — 4.3, the approximate ZAMS 
gravity for stars in the relevant mass range. Stars whose best-fit SED models required this maximum value are indicated 
by log g entries of "4.3", without error bars. The uncertainties listed for stars whose best-fit \ogg values are close to this 
maximum may be somewhat underestimated since the upper range in the error simulations was truncated in the same 
way. 

b The allowed values of v tur b were constrained to lie between and 8 kms -1 , i.e., the range of values available in 
the ATLAS9 model grid. Stars whose best-fit SED models required these limiting values are indicated by v tur b entries 
of "0" or "8", without error bars. The uncertainties for stars with best-fit v tur b values close to these limits may be 
underestimated due to this truncation. 
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Best-Fit Extinction Curve Parameters for Stars in Figures 3-9 
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04 





73 


± 





40 





49 


± 





07 


2 


63 


± 





22 





23 


± 





05 


HD 


142315 


3.37 ± 





34 


1 


28 


± 





21 


1 


90 


± 





18 


1 


32 


± 





02 


-0 


00 


4 


57 


± 





02 





85 


± 





07 





98 


± 





50 





38 


± 





09 


2 


60 


± 





40 





29 


± 





07 


HD 


142378 


3.64 ± 





18 


1 


45 


± 





11 


2 


28 


± 





25 


1 


33 


± 





02 


-0 


00 


4 


64 


± 





02 





03 


± 





08 


1 


51 


± 





40 





41 


± 





08 





91 


± 





22 





09 


± 





06 


HD 


143567 


2.76 ± 





18 





90 


± 





12 


1 


88 


± 





08 


1 


33 


± 





01 


-0 


00 


4 


59 


± 





01 





92 


± 





00 





12 


± 





52 





51 


± 





09 


3 


53 


± 





45 





24 


± 





08 


HD 


145554 


3.65 ± 





19 


1 


40 


± 





12 


1 


79 


± 





10 


1 


31 


± 





01 





00 


4 


58 


± 





01 





95 


± 





05 





76 


± 





40 





33 


± 





07 


3 


54 


± 





34 





24 


± 





06 


HD 


146001 


3.53 ± 





22 


1 


38 


± 





14 


1 


69 


± 





22 


1 


29 


± 





02 


-0 


00 


4 


54 


± 





01 


1 


26 


± 





15 





12 


± 





74 





43 


± 





11 


6 


00 


± 


1 


61 





30 


± 





05 


HD 


146029 


3.56 ± 





27 


1 


40 


± 





17 


2 


11 


± 





17 


1 


34 


± 





01 


-0 


00 


4 


57 


± 





02 





88 


± 





07 


1 


85 


± 





57 





18 


± 





11 


3 


30 


± 





46 





62 


± 





11 


HD 


146416 


2.81 ± 





36 





93 


± 





23 


2 


05 


± 





22 


1 


34 


± 





02 


-0 


01 


4 


53 


± 





02 





89 


± 





08 





45 


± 





82 





44 


± 





16 


3 


08 


± 





57 





76 


± 





15 



a Values of R(V) and km were not determined for IC 4665 stars, i.e., the first group of entries below (HD 161165 through HD 162028). See the discussion in §5.1. For the 
lightly-reddened stars in the third group of entries (HD 142096 through HD 146416), the values of R(V) and km were constrained to follow the relation fc/^ — 0.63i?(V) — 0.84 
from F04. See the discussion in §5.3. 

b The stars HD 37061, HD 37903, HD 147933, and HD 294264 have Johnson R and / photometry available and thus allowed us to solve for the fourth optical spline point O4 at 
7000 A. The values arc: -0.99 ± 0.02, -0.92 ± 0.04, -1.06 ± 0.04, -1.16 it 0.03, respectively. For the rest of the sample, the optical spline is determined only by the points Oi, 
2 , and O3 at 3300 A, 4000 A, and 5530 A, respectively. 



Table 3 

Mean Extinction Parameters for IC 4665 



Parameter 


Weighted 
Mean a 


Sample 
Standard Deviation 13 


RMS of 
Monte Carlo Errors 


Implied 
Intrinsic Scatter d 


Oi 


1.819 


0.219 


0.244 




o 2 


1.306 


0.024 


0.021 


0.012 


o 3 


-0.004 


0.003 


0.001 


0.003 


x 


4.569 


0.015 


0.011 


0.010 


7 


0.922 


0.096 


0.050 


0.081 


Cl 


0.380 


0.475 


0.428 


0.206 


C2 


0.478 


0.108 


0.085 


0.067 


C3 


3.273 


0.736 


0.359 


0.643 


C 4 


0.500 


0.191 


0.052 


0.184 




5.633 


0.709 


0.319 


0.634 



a Mean values computed using 1/cr 2 weighting, with a values as given in Table 2. 

b Standard deviation of the nine measurements shown in Table 2 for each parameter. 

C RMS value of the nine Monte Carlo-based uncertainties listed in Table 2 for each param- 
eter. 

d Computed based on the assumption that the observed scatter (in the 3rd column) is the 
quadratic sum of the measurement errors (in the 4th column) and the intrinsic scatter. 

c Ab ump (= ^) is the area of the Lorentzian-like 2175 A bump for an extinction curve 
normalized by E(B — V). 
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Fig. 1. — Effects of a ±i-class spectral mismatch error on extinction curves derived by the pair method. In 
each of the four groups, the solid curve shows the true shape of a normalized (but arbitrarily offset) extinction 
curve affecting a Bl, B2, B5, or B9 star. The dashed and dash-dot curves show the extinction curves which 
would be derived via the pair method in the presence of a ±| spectral type mismatch error in the cases of 
E{B — V) = 0.30 and E(B — V) = 0.15, respectively. Detailed spectrophotometric extinction curves are 
typically only available in the UV spectral region, which is highlighted in the figure. The extension of the 
mismatch curves into the optical and IR is shown by the dotted curves. Note that the mismatch curves fall 
below the true curve in the UV and above the true curve in the IR when the standard star is cooler than 
the reddened star, and vice versa. The feature in the mismatch curves near 1/A ~ 2.7/im~ 1 is due to the 
Balmer jump. The large feature at 8.2 /im _I (i.e., 1215 A) is due to mismatch of the strong stellar Lyman a 
absorption line. In producing these curves, we use the temperature scale: BO. 5 = 28000 K, Bl = 25000 K, 
B2 = 20000 K, B3 = 18000 K, B5 = 15000 K, B7 = 13000 K, B8 = 12000 K, B9 = 11000 K, and AO = 
9500 K. 
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Fig. 2. — A flexible, parametrized representation of the IR through UV extinction curve. The curve consists 
of two parts: 1) the UV (i.e., A < 2700 A, solid curve in the figure) where we adopt the 3-component 
parametrization scheme of FM90; and 2) the optical/IR (i.e., A > 2700 A, dashed curve) where we adopt 
a cubic spline interpolation through sets of IR (I1J2, I3, optical (0\, O2, O3, O4), and UV (Ui, U2) 
"anchor points." The values of the anchor points and the six parameters describing the UV portion of the 
curve are determined by fitting the observed SED of a reddened star, as described in §3. The dotted line 
shows the "linear component" of the UV extinction curve. The other two components, which are added to 
the linear component, are the 2175 A bump profile centered at A -1 ~ 4.6ttm~ 1 and the far-UV rise, which 
refers to the departure of the curve from the extrapolation of the linear component at very short wavelengths. 
See the discussion in §3.2. The particular curve shown in this figure corresponds to that derived in §5.2 for 
the star HD 147933. 
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Fig. 3. — SED fits for nine lightly-reddened mid-to-late main sequence B stars in the open cluster IC 4665. 
Histogram-style curves are the best-fitting, reddened ATLAS9 stellar atmosphere models, arbitrarily shifted 
vertically for clarity. In the UV (A < 3000 A), small filled circles are IUE low-dispersion spectrophotometry 
binned to match the resolution of the models. In the optical, circles, triangles, and diamonds are Johnson 
UBVRI, Stromgrcn uvb, and Geneva \JB\B-i VV\G data, respectively. In the IR, the large open and closed 
circles are Johnson JHK and 2MASS JHK data, respectively. The photometric data have been converted to 
flux units for display purposes only. As discussed in §5, the comparison between models and photometric 
magnitudes, colors, and indices is performed in the native format of the photometry. 
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Fig. 4. — Extinction curves derived using model atmosphere calculations for nine mid-to-late B main sequence stars 
from the open cluster IC 4665. The values of E(B — V) range from 0.14 to 0.25 for this group. The filled symbols 
show IUE spectrophotometry in the UV (A -1 > 3.33/im _1 ), Johnson, Stromgrcn, and Geneva photometry in the 
optical, and 2MASS JHK in the near-IR. Open symbols indicate near-IR Johnson JHK data. The solid curves are 
the fits to the data (i.e., the "flexible, parametrized curves" described in §3.2) as determined by the SED-fitting 
procedure, and the shaded regions show the 1-cr uncertainty in the curves, based on Monte Carlo simulations. For 
comparison, the dash-dot curves show the average Galactic extinction curve (corresponding to R(V) = 3.1). Because 
of an apparent conflict in the 2MASS and Johnson JHK data, we assumed that the shape of the IC 4665 curves in 
the near-IR follows the average Galactic curve. The top solid curve and the square symbols show the simple mean 
of the nine IC 4665 extinction curves, along with the average Galactic curve for comparison. The error bars indicate 
the sample standard deviation of the individual curves. Note that this scatter is comparable to the 1-cr errors of 
the individual curves, indicating that the shape of the extinction curve in IC 4665 is uniform to within our (small) 
measurement errors. 
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Fig. 5. — Comparison between extinction curves derived for the cluster IC 4665 using the extinction- without- 
standards technique described in this paper (upper set of curves) and using the traditional pair method 
approach (lower set of curves; from HHT). The average Galactic curve is shown with both sets of curves for 
reference (dash-dot curves). 
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Fig. 6. — Same as Figure 3, but for a group of nine moderately reddened early B stars. 
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Fig. 7. — Same as Figure 4, but for nine moderately reddened early B stars. Dotted curves show pair method 
extinction curves for all stars except HD 294264. Those for HD 210121 and HD 27778 are unpublished results 
by us, and the rest are from the catalog of FM90. We are not aware of any pair method extinction curves 
for HD 294264. The values of E(B — V) range from 0.36 to 1.08 for this group. These curves illustrate the 
wide range observed in the UV-through-IR wavelength dependence of interstellar extinction. 
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Fig. 8. — Same as Figure 3, but for a group of nine lightly-reddened mid-to-late B stars in the direction 
347° < I < 355° and 18° < b < 26°. 
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Fig. 9. — Same as Figure 4, but for a group of nine lightly-reddened mid-to-late B stars in the direction 
347° < I < 355° and 18° < b < 26°. The values of E(B - V) range from 0.09 to 0.22 for this group. The 
dotted curve, superimposed on the average curve at the top of the figure shows the result for the nearby star 
HD 147933, located at (I, b) = (353.7°, 17.7°) with E(B - V) = 0.50 (see Figure 7). 
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